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Abstract 

We present a calculation of the finite volume corrections to meson masses and 
decay constants in three flavour Partially Quenched Chiral Perturbation Theory 
(PQChPT) through two-loop order in the chiral expansion for the flavour-charged 
(or off-diagonal) pseudoscalar mesons. The analytical results are obtained for three 
sea quark flavours with one, two or three different masses. We reproduce the known 
infinite volume results and the finite volume results in the unquenched case. The 
calculation has been performed using the supersymmetric formulation of PQChPT 
as well as with a quark-flow technique. 

Partial analytical results can be found in the appendices. Some examples of cases 
relevant to lattice QCD are studied numerically. Numerical programs for all results 
are available as part of the CHIRON package. 



1 Introduction 


Quantum chromodynamics (QCD) is nowadays accepted to be the theory describing the 
strong force. The smallness of the coupling constant at high energies makes it possible to 
test and conhrm the theory in highly energetic scattering. It also provides - at least in 
principle - a way to obtain various low-energy hadronic observables, such as masses and 
decay constants, but it has hitherto been impossible to derive such quantities of interest in 
terms of analytical expressions by means of ab initio calculations. A numerical approach 
that can circumvent the problem is lattice QCD. A review of the applications to flavour and 
low-energy hadron physics is [1] . To calculate observables, one uses a numerical evaluation 
of the QCD path integral in a Monte Carlo approach. A number of restrictions follow 
from the nature of the calculation. Since it is carried out on a space-time lattice in a 
hnite volume, it is of high interest to have the effect of the hnite volume under good 
control. Furthermore, although lattice computations in meson physics are now feasible 
when using physical parameters for the light quark masses, a lot of calculations still use 
unphysically high masses for the quarks. It is also useful to vary quark masses to study a 
number of phenomena. A common solution to study quark mass dependence with lower 
computational needs is given by partial quenching. In partially quenched QCD (PQQCD), 
one associates different masses (usually larger ones) to the sea quarks and the valence 
quarks. Valence quarks are those connected to the external operators while sea quarks 
are those in the fermion determinant or equivalently in closed loops. Sea quarks are only 
connected to external states via gluons. 

The preferable way to correct for unphysical quark masses is by means of Chiral Per¬ 
turbation Theory (ChPT) |2],|3lll]. Finite volume effects for ChPT have been introduced 
in mum- The corresponding effective theory for PQQCD is given by Partially Quenched 
Chiral Perturbation Theory (PQChPT) [8]. The arguments underlying this are elaborated 

in [9]. 

The proper matching of calculations in PQChPT to results from Partially Quenched 
Lattice QCD allows a whole new landscape of possibilities, such as improved validation and 
extrapolation of lattice results, or a more accurate determination of the chiral low-energy 
constants (LECs), see e.g. [10]. It should be stressed that, as opposed to fully quenched cal¬ 
culations, partially quenched calculations are connected to their corresponding unquenched 
scenarios by a continuous change in variables, making it possible to immediately extract 
physical results from otherwise unphysical simulations. 

In this paper, we address the hnite volume corrections through two-loop order in the 
PQChPT framework, specihcally for the Favour-charged or off-diagonal mesons. The in- 
hnite volume (IV) results in PQChPT to this order are known for three [TTl fl^ [T3] and 
two [13] sea quark Favours. The Fnite volume (FV) corrections in (unquenched) ChPT at 
two-loop order have been addressed in our earlier study [13]. The needed integrals have 
been worked out in [TB]. Our expressions are valid in the frame with vanishing spatial 
momentum, p = 0, often called the center-of-mass frame. In the so-called moving frames 
or with twisted boundary conditions there will be additional terms. We have chosen to 
present our result in terms of lowest order masses given the ambiguity in expressing the 
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results in terms of the large number of possible different physical masses. 

Earlier work on hnite volume corrections at NNLO are besides our own work [15], the 
pion mass in two-flavour ChPT im and the vacuum expectation value in three flavour 
ChPT [TS|. Extensions of the latter work to partially quenched are in [12]. We did not 
hnd published results for the hnite volume corrections at one-loop order in the partially 
quenched case. They are however implicit in the expressions given for the staggered par¬ 
tially quenched case in [20] [21] . 

We give a short list of references for ChPT and discuss some small points in Sect. The 
dehnitions of the integrals we use and how they relate to the results in [TB] is given in Sect. 
The next section describes our major result which is the full hnite volume correction to the 
pion mass and decay constant to two-loop order in ChPT. Sect, [^contains the results for 
the three-havour case for pion, kaon and eta for both the mass and decay constant. The 
large two-loop order formulas are collected for one case in the appendices and all of them 
can be downloaded from [22] • A numerical discussion of our results is in Sect. 


2 Partially Quenched Chiral Perturbation Theory 

This section is very similar to the description of PQChPT given in [TB] since our work is 
the extension to hnite volume of that paper. 

An introduction to ChPT can be found in [23] and in the two-loop review [23] • The 
lowest order and p^-Lagrangian can be found in [3]. The order p® Lagrangian is given in 
[25] . We use the standard renormalization scheme in ChPT. An extensive discussion of the 
renormalization scheme can be found in [26] and Important for our work is that the 
LECs do not depend on the volume [7]. An introduction with applications to lattice QCD 
is [2E]- References to more introductory literature can be found on m 

The expansion in ChPT is in momenta p and quark-masses. We count the latter as two 
powers of p. This counting is referred to as p-counting. We prefer to designate orders by 
the p-counting order at which the diagram appears. Thus we refer to lowest order (LO) as 
order p^, next-to-leading order (NLO) as order p^ or one-loop order and next-to-next-to- 
leading order (NNLO) as order p® or two-loop order and include in the terminology one- 
or two-loop order also the diagrams with fewer loops but the same order in p-counting. 


2.1 The Lagrangian 

Three massless quark flavours QCD has a chiral symmetry 

G = SU{nf)LxSU{nf)n. ( 1 ) 

which is spontaneously broken to the diagonal subgroup SU{3)v- The Goldstone bosons 
following from this spontaneous breakdown are described by the meson octet matrix 


/ 


0(a;) = 


P2 


TT 

K- 




K+ 






( 2 ) 


■ 73 ^ 
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The flavour-singlet component has been integrated out since it is heavy due to the U{1)a 
anomaly. The spontaneous symmetry breaking is the basis of ChPT. 

In partially quenched QCD one distinguishes between valence and sea quarks. Valence 
quarks are connected to the external states (or operators) while the sea-quarks are those 
contributing in closed loops only connected via gluons to external states. These can be 
given different masses in lattice QCD calculations. The ChPT for this partial quenching 
can be done by studying the quark flow generalizing the quenched case studied in |29j. 
One can then treat the the sea and valence lines differently. Alternatively, one can make 
use of the supersymmetric formulation of PQChPT [S]. In the latter, three corresponding 
sets of quarks are introduced instead of only two: In addition to the valence and sea sector, 
a set of so-called ghost quarks is added. These are “bosonic” in the sense that they are 
treated as commuting variables. With their masses fixed to the same numerical values 
as present in the valence sector, they will cancel exactly the contribution coming from 
closed valence quark loops. Most of the remainder of this section will be concerned with 
the supersymmetric formulation. The changes needed to use a quark flow technique are 
discussed at the end. 

The chiral symmetry group is formally extended to the gradecQ 

G = SU (rival + ^sea|nval)L X SU (rival + nsea|^val)il , (3) 

for the case of rival valence and risea quarks. The chiral group G is spontaneously broken 
to the diagonal subgroup S'!/(rival + ’^sea|’^vai)v- We will work in the flavour basis rather 
than in the meson basis. We will thus use fields (pat corresponding to the flavour content 
of qaqb- The mixing of the neutral eigenstates and the integrating out of the singlet degree 
of freedom is taken care of by using a more complicated propagatoi0 

The corresponding Goldstone degrees of freedom are put in a matrix with the generic 
structure 



/ 

qvQv 


qvqs 


qvqB 

\ 


$ = 


qsQv 


qsqs 


qsqB 


(4) 


1 

qBqv 


qsqs 


qBqB 

) 



V denotes valence, S denotes sea and B denotes the bosonic ghost quarks. Note that the 
meson fields containing one single ghost quark only will themselves obey fermionic, i. e. 
anticommuting, statistics. 

The structure of the Lagrangian is similar to standard ChPT for a generic number of 
flavours. The lowest order Lagrangian is 

A = ^(V^ + X+)- (5) 

Mhe precise structure of the symmetry group is somewhat different, but the one given here is sufficient 
for both the present discussion as well as for practical calculations. The “approximate” symmetry group 
reproduces the right Ward identities [iniiso]. 

^This is described in detail in m- It is possible to use the same method also in standard ChPT. 
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At one-loop, it is given by 

£4 = Lo {U^U^U^UJ) + Li + 1/2 {uV){Uf,U^) + Z/3 

+1/4 {u^Uy){x+) + -^5 + Lq (x+)^ + i /7 (x-)^ + {x+ + xi) + • • • (6) 

We show only the terms relevant for our work. 

The generalized Goldstone manifold is parametrized by 

u = exp (i^ /{\/2F)'j ( 7 ) 

similar to the exponential representation in standard ChPT. It is a 9 x 9 matrix with 
fermionic parts. We have furthermore introduced 

X± = u^xu^ u. (8) 

The matrix y is for this work restricted to 

X = 25odiag(mi,... , 7719 ) (9) 

with rrii the quark mass of quark i and Bo a LEG. We have here mi = mj, m 2 = mg, m 3 = 
mg as the valence masses and m 4 ,m 5 ,m 6 as the sea quark masses. Ordinary traces have 
been replaced by supertraces, denoted by ( ), defined in terms of the ordinary ones by 

Str(^ ^)=TrA-TrT). (10) 

B and C denote the fermionic blocks in the matrix. The supersinglet $ 0 ) generalizing the 
77 ', is integrated out to account for the axial anomaly as in standard GhPT, implying the 
additional condition 

($)= Str (<!))= 0 . (11) 

However, as mentioned above, we will work in the flavour basis enforcing the constraint 
0 via the propagator. 

A calculation in PQGhPT has to be performed using a larger set of operators since no 
further reduction by means of Gayley-Hamilton relations can be performed. The three- 
flavour PQGhPT Lagrangian (equation (|^) thus has 11 LEGs for PQGhPT. 

The LEGs for standard three flavour GhPT are related to those of three flavour PQGhPT 
via 

Ll = L{ + Ll/2, LI = Ll + LI, L; = 1^-2 L^, (12) 

and LI = for the others. Note that a numerical value for Lg cannot be obtained by 
experiment, but can be determined only via PQQGD lattice simulations or modelling. 

An additional comment is that the divergences for PQGhPT are directly related to 
those for risea-flavour GhPT IZ71 when all traces are replaced by supertraces. This can be 
argued using the formal equivalence of the equations of motion used or via the replica trick 

Ell. 
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2.2 The propagator and notation for masses and residues 


The variant of PQChPT, considered in this paper, conies with three valence qnarks, with 
masses mi, m 2 , m 3 and three sea qnarks with masses The additional ghost 

qnarks emerging only in the snpersymmetric formnlation have masses m 7 ,m 8 ,m 9 . They 
do not appear explicitly since they are hxed to the ones in the valence sector, i.e. my = mi, 
ms = m2, mg = m3. 

We nse the nnmbers d^ai and dgea to denote the nnmber of non-degenerate qnark masses 
in each sector. In the case of two non-degenerate mass scales for one sector, it is the 
two masses with the lowest indices that we set degenerate, which will in tnrn both be 
represented by the mass scale with the lowest index, e.g. in the case dsea = 2 we have 
mi = ms 7 ^ me and expressions will be explicitly dependent on mi and me only. 

In fact we will always absorb a factor 2Bq in the notation and we nse 


Xi = 2Bami, xtj = i,(Xi + Xi) ■ 


(13) 


The lowest order masses for off-diagonal mesons with flavonr content qiq^ are given by 
Xij and we will use Xi rather then xa for equal masses. Dealing with masses for the 
diagonal valence mesons in PQChPT is not trivial. This is discussed in detail in |30] and 
extended to NNLO in [32]. The diagonal sea quark sector has two masses associated with 
it, corresponding to the neutral pion and eta masses. These we denote by Xtt and Xv They 
are dehned as the solutions to the equations 


Xtt T Xr) 


X'lrX'n 


X (X4 + X5 + Xe) ; 
^(X4X5 + X5Xe + X4X6). 


(14) 


They are non-polynomial in the sea masses Xj for three non-degenerate quark masses, i.e. 
dsea = 3. For dsea = 2 oue has instead = X 4 and Xv = (1/3) (X 4 + 2x6)- 

The flavour-charged propagator, connecting with 0^,, is given by [H HOl [30| 

= ,2 '' 4, . (15) 

- Xij + 


with Xij = iXi + Xi)/ 2 , the lowest order meson mass, and the signature ej is dehned as -t -1 
for the havor indices of the + risea fermionic quarks, and as —1 for the havor indices of 
the rival bosonic ghost quarks. In the present calculation, with the number of valence and 
sea quarks as given above, ej thus takes the values 


e 


j 


-f 1 for j = 1 ,..., 6 
-1 for j = 7,8,9. 


(16) 


The havour-neutral propagator, connecting a havour held 0** to (jyjj, on the other hand 
suhers from additional contributions emerging from the elimination of the $0 and the 
partial quenching [HI [10], [30] . We write it as 


Gl(k) = G%(k)S„-G%{k)ln. 


(17) 
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The additional terms are either 


iGUk) = 


K 


’JTTT) 


+ 


dJ 


~ Xi + k"^ ~ Xj + 


+ 


K 


m3 


jdV 

■^Tvij 


k"^ - Xn + i£ k"^ - Xri + i^' 


for the case with i ^ j and Xi 7^ Xj) or 

-iGlAk) = 


Rf 


+ 


Rf 


{k'^ — Xi + k‘^ — XiR R 


R. 


r]ii 


— + 


tdV 

■^nii 


k"^ - XiT + k? - Xri + iA 


(18) 


(19) 


for the case with Xi = Xj which clearly includes i = j. In the second case, the sum of single 
poles is supplemented with an unphysical double pole. Since double poles emerge due 
to the partial quenching in the valence sector, they disappear by taking the appropriate 
unquenched limit. 

Using the ratios of products of differences of masses 


r:^ = 

Xa Xb) 


Rabc = 

Xa Xb 


Xa Xc 


TDZ _ 

^abcd 

iXa Xb){Xa Xc) 


Xa-Xd 


)Z 

(Xa Xb){Xa Xc){Xa 

- Xd) 

ahcdefg 

iXa XejiXa Xf){Xa 

-XgV 


( 20 ) 


the residues R of the neutral meson propagator in equations (18) and (19) are (for dgea = 3) 

( 21 ) 


TDl _ JDZ 

^jkl 


ryd _ -QZ 

— -^24567r7?5 


7DC _ r^i I Tyi _j_ -ni _ T^i _ T^i 

-^57777 ' ^TrrjT} ^TTirr}' 


Note that many of these quantities vanish when i takes the value of a sea quark index. 
The sea-quark propagators thus do not contribute any double poles as expected since these 
originate from the quenching in the valence sector. 

For dsea = 2 or = xs = X 4 . The needed residues simplify to 


p* _ pz 

^jk 


prf _ pz 


pc _ pi \ pi 

nij — 44 , 4 ^ -TLq^ 


R 


rjT) 


( 22 ) 


The corresponding propagator can be obtained by removing all pion indices as well as the 
pion mass pole from equations (18) and (19). 

The physically less interesting case dsea = 1 immediately yields Xn = Xv — Xe = Xb = 
X 4 - All residues from the sea quark sector are reduced to numbers, only 


pi _ pz pd _ p 

rty — n^4j, rtj — rt. 


i4 ) 


(23) 


appear. 
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2.3 The quark flow case 

We have performed the calculation using the supersymmetric method described above but 
also with the quark flow method [22]. We use the same Lagrangians as in ([^ and (|^ but 
with normal traces everywhere. The matrix $ is now written in terms of generic helds (pij 
and all indices are kept symbolic implying summations. 

Connecting propagators of a held 0*^ to should be done by using 

Gijki{k) = G^j{k)6ii6jk — 5ij5kiG%{k)/Usea. ■ (24) 

The propagators Gij{k),Gjj{k) remain the same but we can now disregard the factors e^ 
since with this method there are no bosonic ghost quarks. 

After constructing the Feynman diagrams using the above, the quark how is visible 
following the symbolic havour indices. Next, one replaces the index lines that connect to 
external helds or operators by their appropriate valence value. The remaining index lines 
are now sea indices and are summed over with the sea quark indices. 

The results obtained with the quark how method agreed in all cases with those of the 
supersymmetric method. 


3 The finite volume integrals 


The loop integrals at hnite volume at one-loop are well known. There is a sum over 
discrete momenta in every direction with a hnite size rather than a continuous integral. 
The Poisson summation formula allows to identify the inhnite volume part and the hnite 
volume corrections. The remainder can be done with two diherent methods. For one- 
loop tadpole integrals the hrst method was introduced by mum and a sum over Bessel 
functions, that for large ML converges fast, remains to be done. With the other method 
one remains instead with an integral over a Jacobi theta function, this method can be used 
for small and medium ML as well. It can be found in [33]. The extensions to other one- 
loop integrals is done in both cases by combining propagators with Feynman parameters. 
The hrst method was extended to the equal mass two-loop sunset integral ini and later to 
the more general mass case in [15]. The latter extended the Jacobi theta function method 
as well to the sunset case. Details and further references can be found in [TB]. In this paper 
we use Minkowski notation for the integrals. 

For the one-loop integrals needed here, we use a notation that does a hrst classi- 
hcation according to the sum of the powers of the propagators with diherent masses, 
mi, m 2 ,..., mrnaa:- We label the integrals A,B,C,D for a total power of propagators of 
n = 1,2, 3,4 respectively, since total powers of up to 4 can appear in the calculation as 


follows from the discussion of double poles in Sec. 2^ The diherent mass scales are given 
as consecutive arguments of the integral. Alternatively, if only one mass scale in total is 
present, we omit its repetition as a shorthand notation. For the present calculation at most 
two diherent scales can appear. 
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Both scalar and tensor integrals will occnr, e. g. in the simplest case of one single 
propagator raised to single power 




d‘^r 


'v {2 'kY 




(25) 


We nsed the snbscript V to indicate it is a hnite volnme snm and integral. 

More Lorentz strnctnres are possible than in the inhnite volnme case. We dehne the 
tensor as the spatial part of the Minkowski metric to express these. For the 
center-of-mass (cms) case this is snfficient. The needed fnnctions for the above example 
are 

= gt,uA22{rrY) + • (26) 

We then nse Passarino-Veltman identities in order to fnrther simplify the resnlt. In inhnite 
volnme the relation obtained by considering can be nsed to remove A 22 - In 

hnite volnme, we again remove the 2422 -type integrals from the extended relation 


dA 22 {m‘^) -l- 3A23{m‘^) = m?A{m?). 


(27) 


Each integral is split into an inhnite volnme contribntion and a hnite volnme correction 
by means of the Poisson snmmation formula, while simultaneously being expanded in e up 
to the necessary order. 


2 

A{nY) = ^0^^ + + e (^AYnY) + A'^YmY)) H-• (28) 

Here, Aq = ^ -l- log( 47 r) + 1 — 7 . The same split is done for all one-loop integrals. The 
expressions can be obtained by using the relations 


B{nY) 


C{rrY) 


D{nY) 



A{ml) — 24 ( 7712 ) 

9 91 

mf — 7772 

The sunset integrals, dehned as 


( U LJ LJS TT \ 

fliy) fll/ j 

{n, ml, ml, ml, p) = 


Y r d'^r d'^s _ {1, r^, r^r^, _ 

Jv { 27 rY { 27 rY (r^ — mf)”^ (s^ — ((r -f- s — pY — ’ 


(29) 


(30) 


now come with eight diherent pole conhgurations. We label these by the index n according 
to Tab. [^analoguous to the inhnite volume dehnitions of [HI [121 [131 IE]- 










77i 772 ns 

77 = 1 

111 

77 = 2 

2 11 

CO 

12 1 

(77 = 4) 

112 

77 = 5 

2 2 1 

(77 = 6) 

2 12 

77 = 7 

12 2 

00 

2 2 2 


Table 1: Overview of the notation for the possible conhgurations of powers of propagators 
in the H functions in PQChPT. Redundant conhgurations are given in parentheses. 


The interchange (r, O ( 5 , 7712 ,^, 2 ) allows to show that are related di¬ 

rectly to if)), if". can also be related to if^^ using the trick shown in [M] and also 
used in [16], now taking the pole conhgurations into account properly. The resulting if^i. 
and if^ can then be reduced to six pole conhgurations only, cf. table the bracketed 
ones can be eliminated via the interchange above. In the scalar case if, only four pole 
conhgurations are needed. 

For the partially quenched calculation we thus generalized the sunset integrals used in 
our earlier work via 

H{Xi,Xj,Xk]p'^) ^ H{n,Xi,Xj,Xk]p‘^), (31) 

introducing the new index n for the pole conhgurations as the hrst argument. Note on the 
side that all new pole conhgurations are related to the simplest one by diherentiation with 
respect to the mass scales. 

In the cms frame, we reduce the tensor structure of the sunsets as 


if^ = p^Hi (32) 

Hfiv = P^PuH21 + gfiuH22 + t^i/if27 • 


As in [T5|, we renormalize the FV sunsets by not only subtracting the inhnite part but also 
an additional hnite part containing 0{e) contributions of one-loop integrals. In this way, 
the latter integrals will cancel out of the hnal result, and thus do not need to be computed. 
The splitting for 77 = 1 


m = 


An 




(ml) + (ml) + A^(ml)^ + (^A^^{ml) + A^^{ml) + A 


Ve, 


7T7c 
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Figure 1: Diagrammatic contributions to the pseudoscalar self-energy, up to 0{p^). Cir¬ 
cular vertices are of 0{p‘^), the hlled boxes are of 0{p‘^), the open box is of 0{p^). The 
tree level diagrams (a,b,i) do not contribute to hnite volume corrections. 



has to be generalized for the other pole conhgurations by taking the appropriate derivatives 
w.r.t. the masses. 


4 Analytical results 

The calculation of the masses proceeds in the usual way from the Feynman diagrams for 
the self-energy shown in Fig. [1} We have performed the calculation for the off-diagonal 
mesons, i.e. consistening of a valence quark and a different valence anti-quark, and for the 
case of three flavours of sea quarks. The calculation has been done for all mass cases, equal 
and different valence quark-masses, d^ai = 1, 2, and sea quark masses all equal, dgea = 1, 
two equal and the third different, dsea = 2 and all three different, dgea = 3. 

A large number of checks have been done on the calculations. They have been performed 
both in the supersymmetric formalism and using quark flow techniques. The inhnite volume 
results are also in full agreement with [mini [ 13 ]. The hnite volume parts agree with our 
earlier results na when these are expressed in terms of lowest order masses and when the 
sea masses are put equal to the valence masses. 
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The formulas especially for the case of three different sea quark masses are very long. In 
App. I^we list the case of equal valence masses and two sea quark masses. This corresponds 
to the charged pion mass in the isospin limit. The other cases can be downloaded from 


The masses are given as 


2 I 2(4) 

^ij = Xij + "iq 


+ A 


V 2(4) , 2(6) 

rrii] + rrii] 


+ A 


V 2(6) 

m 




(34) 


In addition a superscript indicating d^aidgea is added. The inhnite volume and the one- 
loop hnite volume corrections were known before. The new parts are the two-loop hnite 
volume corrections. These we split in addition in an dependent part and a pure two-loop 
contribution 

AUnf > = + A*-. (35) 

The subscript ij is set to 12 for d^ai = 1 and to 13 for d^ai = 2 similar to the inhnite volume 
work. 

The decay constant is dehned in the usual way as 

= iV2FijPf ,, (36) 


for the pseudoscalar meson Mij with quark content i ^ j and momentum p. The calculation 
needs the the diagrams of Fig. [^for the wave function renormalization and the same ones 
with one external meson leg replaced by an insertion of the axial current. 

We split the result as 

Fij = Fo + F,f + A^F^f + fIj'> + A^F,f . (37) 

The NNLO part is split again in 

A^F,f = A^F^-^^ + A^Ffl ^^^. (38) 

The calculations have been done using the supersymmetric and the quark how methods. 
The inhnite volume and NLO results agree with the known expressions and the result 
reduces in the correct limit to the unquenched results of our earlier work [15] . The formulas 
are rather long, the case for equal valence masses and two diherent sea masses corresponding 
to the charged pion decay constant in the isospin limit is given in App. The expressions 
for the other cases can be downloaded from izg. 


5 Numerical examples 

The intention is that various lattice QCD collaborations can use our formulas. All cases 
discussed have been included in the package CHIRON [35] available from [36]. The nu¬ 
merical results shown in this section have been obtained with that implementation. The 
programs have been cross-checked with an independent version. It has been checked that 
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the results reduce in the appropriate limits to those of our earlier work [15]. For this pur¬ 
pose the expressions obtained in [15], but rewritten in terms of lowest order masses and 
decay constants, have been implemented and included in CHIRON [3S]- In addition, a 
check has been done that the different mass cases reduce to each other numerically. 

For input values we have chosen the recent global £t for the n EH. We have set the 
extra LEG Lg = 0. We always use a scale of /i = 0.77 GeV. For the size of the lattice we 
present results for a length L such that ML = 2 for M = 0.13 GeV. The lowest order pion 
decay constant we have chosen throughout as Fq = 87.7 MeV. 

The numerical results are presented via 




m 


2V 

ij 


m 


2oo 


Xij 

pV _ poo 

P n 

Fo 


(39) 


We thus plot the size of the hnite volume corrections relative to the lowest order value 
of the quantity under consideration. Note that the results are for charged or off-diagonal 
mesons. They consist of a quark and a different anti-quark which might have equal mass. 


5.1 dval — ^sea — 1 

Here we set all valence and all sea masses equal, d^ai = dgea = 1- The size of the hnite 
volume corrections as a function of xi and X 4 is shown in Fig. [^ The corrections in this 
case are reasonable, at most a few %, except for very low masses and become very large 
for low valence and high sea quark mass. 


5.2 The pion mass 

In this subsection we look at the case where the lowest order mass is around the pion 
mass. We plot with y/Xi 2 = 0.13 GeV. The strange sea quark mass we have always 
chosen such that the average lowest order kaon mass is 0.45 GeV. This corresponds to 
yjx^ = ^2(0.45)^ — (0.13)^ GeV ~ 0.623 GeV. The other input parameters are chosen as 
given in the introduction of this section. We have restricted the sea up and down quark 
masses corresponding to a lowest order sea quark pion of 100 to 300 MeV. 

The hrst case we look at is dvai = 1, dsea = 2. This corresponds to taking the up and 
down quark masses equal in both the valence and sea quark sector and a different strange 
quark mass. This is the isospin limit. The result is shown in Fig. |^a). There is a rather 
large cancellation between the and p® correction while the p® contribution coming from 
the L\ is fairly small. 

We now include isospin breaking in the valence sector. We thus look at the case 
with dvai = 2, dgea = 2. We £x the valence quark masses such that Xi + X 2 = 2xi2 and 
Xi/X 2 = 1/2. There is a sizable isospin breaking visible in the hnite volume corrections, 
as shown in Fig. [^b). 


12 






Figure 2: The finite volume corrections relative to the lowest order value as defined in 
(39) for the case with all valence masses equal and all sea masses equal. Left: the 

correction to the mass-squared, contour lines are drawn at 0.03,0.01, 0.003, 0, —0.03, —0.01 
starting from the bottom left and going counterclockwise. Right: The correction to the 
decay constant, contour lines are drawn at —0.001, —0.002, —0.005, —0.01, —0.02, —0.05 
going from top-right to bottom-left. 


The opposite case, isospin breaking in the sea sector, but not in the valence sector, leads 
to numerically similar but opposite sign corrections. Here we used Xi = X2, Xi = X 5/2 and 
X 4 -f X 5 = 2xav The results are shown in Fig. [^c). 

Finally, we introduce isospin breaking in both the valence and sea quark sector with 
Xi/X 2 = 1/2, Xi = X 5/2 and Xi + X 5 = 2Xav The results are shown in Fig. |^d). The total 
isospin corrections are rather small. 

The numerical cancellation between the isospin breaking in the valence and sea quark 
case is accidental. The corrections due to valence and sea quark masses are all second order 
in isopin breaking. The same argument as in the unquenched case goes through both for 
the valence and sea quark masses. We have compared four scenarios in Fig. We show 
the and the full -|- p® result first with no isospin breaking, then only in the valence 
sector or only in the sea sector and finally in both sectors. The curves are those shown in 
Fig-i a-d). We have checked numerically by using a different ratio for the isospin breaking 
that the corrections are indeed second order in isospin breaking. 

5.3 The pion decay constant 

In this subsection we look at the same cases as before. The lowest order mass is around the 
pion mass. We plot with yTi 2 = 0.13 GeV. and as before yTe = ^2(0.45)^ — (0.13)^ GeV 
~ 0.623 GeV. The other input parameters are again chosen as given in introduction of this 
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(c) (d) 

Figure 3: The corrections for the pion mass relative to the lowest order mass as a function 
of the average up and down sea quark mass via Xav (a) The isospin limit, Xi = X 2 , 
Xi = X 5 = Xav (b) Isospin breaking in the valence sector, xi = X 3/2 and X 4 = X 5 = Xav 
(c) Isospin breaking in the sea sector, xi = X 2 and X 4 = X5/2- (d) Isospin breaking in both 
sectors, xi = X 3/2 and Xi = X 5 / 2 . 
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Xav [GeV^] 


Figure 4: Comparing the finite volume correction for the meson masses for the cases with 
no isospin breaking (none), only in the valence sector (val), only in the sea sector (sea) 
and in both (full) for the meson mass squared. The upper curves are the the bottom 
the + p® results. 


section. We have restricted the sea up and down quark masses corresponding to a lowest 
order sea quark pion of 100 to 300 MeV. 

The first case we look at is dvai = 1, <^sea = 2. This corresponds to taking the up and 
down quark masses equal in both the valence and sea quark sector and a different strange 
quark mass, i.e. the isospin limit. The result is shown in Fig. [^a). The total p® correction 
is fairly small. 

We now include isospin breaking in the valence sector. We thus look at the case 
with dvai = 2, dgea = 2. We fix the valence quark masses such that Xi + X 2 = 2xi2 and 
Xi/X 2 = 1/2. There is a sizable isospin breaking visible in the finite volume corrections, 
as shown in Fig. [^b). 

The opposite case, isospin breaking in the sea sector but not in the valence sector leads 
to numerically much smaller effects. Here we used Xi = X 2 , X4 = X5/2 and Xi + Xb = 2 Xav 
The results are shown in Fig. |^c). 

Finally, we introduce isospin breaking in both the valence and sea quark sector with 
Xi/X 2 = 1/2, Xi = X 5/2 and Xi + X 5 = 2Xav The results are shown in Fig. |^d). The total 
isospin corrections are failry small. 

The corrections due to valence and sea quark masses are all second order in isospin 
breaking. The same argument as in the unquenched case goes through both for the valence 
and sea quark masses. We compare the same four scenarios as for the pion mass, no isospin 
breaking, only in the valence sector, only in the sea sector and in both sectors. The curves 
are those shown in Fig. [^a-d). 
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Xav [GeV^] Xav [GeV"] 


(c) (d) 

Figure 5: The corrections for the pion decay constant relative to its lowest order value 
as a function of the average up and down sea quark mass via Xw (a) The isospin limit, 
Xi = X 2 , X 4 = X 5 = Xav (b) Isospin breaking in the valence sector, xi = X 3/2 and 
X 4 = X 5 = Xav (c) Isospin breaking in the sea sector, xi = X 2 and X 4 = X5/2- (d) Isospin 
breaking in both sectors, xi = X 3/2 and X 4 = X 5 / 2 - 
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Xav [GeV^] Xav [GeV^] 

(a) (b) 

Figure 6: Comparing the finite volume correction for the meson decay constant and masses 
for the cases with no isospin breaking (none), only in the valence sector (val), only in the 
sea sector (sea) and in both (full) for the meson mass squared, (a) (b) +p®. 


5.4 The kaon mass and decay constant 

We now look only at the dvai = dsea = 2 case but choose the valence masses such that we 
have a lowest order pion mass of 130 MeV and a lowest order kaon mass of 450 MeV. This 
corresponds to yTi = 130 MeV and yTa ~ 623 MeV. We plot the finite volume corrections 
relative to the lowest order value of the quantity in Fig. as a function of X 4 = X 5 - For 
the sea quark strange mass we use Xe = 1-02x3. The LECs are again the ones from [37] 
and L such that ML = 2 for M =130 MeV. 

For the kaon we see that we reproduce the results of [15] that near the physical case 
the p^ corrections are very small. The total finite volume corrections to the mass remain 
fairly small. The kaon decay constant has larger corrections but they remain in the few % 
region for the parameters considered. 


6 Conclusions 

We have computed the NNLO expressions for the masses and decay constants in three- 
flavour partially quenched ChPT for all possible mass cases. The calculation has been 
performed using two different formalisms, quark-flow and the supersymmetric method. 
The known infinite volume expressions have been reproduced. We quoted the expressions 
for the equal valence and two different sea quark masses in the appendices. The other cases 
can be obtained from [22] • 

The numerical work shows finite volume corrections of a similar size as those in the 
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X4 [GeV^] X4 [GeV^] 

(a) (b) 

Figure 7: The finite volume corrections for a valence mass close to the kaon mass relative 
to the lowest order value, (a) the kaon mass squared, (b) the kaon decay constant. 

unquenched case [15]. We have presented some representative numerics. The numerical 
work has been done using C+ + . The programs are available together with the infinite 
volume results in [36]. The analytical work relied heavily on FORM |38j . 
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A Expressions for the mass 

= +R^(xi) (-1/3x1 (2/3xiR^6i')+5''(Xi) (-1/3 Xi Rue.) 

(40) 


= +R^(xi) ( - 64/3 Llxl RU. “ 32/3 L^^Xi Xe Ru. - 64/3 t,Xi X4 R^e. 
+32 T 7 X 1 + 32 TgXi + 8/3 T 5 X 1 Rue. ~ 8/9 LgXi Xe Ru.^ ~ 16/9 LgXi X 4 Rie.^ 

+40/3 Llxl RU. - 8 C^3Xi Rfie. “ 8 L^xl RU. + 20 L^xl + 8 L[xl - 8 L^Xi R?4e. 

“8 TqXi Rue. ~ 16/3 (xe + 2 X 4 ) CgXi Rue. ~ 16/9 (xe + 2 X 4 ) T 4 X 1 Rie.^ 

—8/9 (xe + 2 X 4 ) C 4 X 1 Ru.^ + 8 (xe + 2 X 4 ) -^ 4 X 1 Rue.) + ^^(Xi 4 ) (32 (x 4 + Xi) -^sXi 


18 














-16 (x4 + Xi) Llxi + 20 (X4 + Xi) Llxi + 8 (x4 + Xi) ^Sxi) + ^^(Xie) (l6 (xe + Xi) ^sXi 

“8 (xe + Xi) -^5X1 + 10 (xe + Xi) -^3X1 + 4 (xe + Xi) -^0X1) + ^^(X4) (48 L^xi Xi 

-48 Llxi Xi + 12 ^ 2 X 1 X4 + 48 L[xi Xi) + ^^(X46) (32 (xe + X 4 ) ^eXi - 32 (xe + X 4 ) ^ 4 X 1 

+8 (xe + X 4 ) -^ 2 X 1 + 32 (xe + X 4 ) -^ 1 X 1 ) + iXv) (l28/3 -^gXi 

+ 32/9 L5X1 Xe '^14»7 ^rj&l 32/9 L^Xi X4 -^^61 + 16 -^ 3 X 1 Xr) -R^ei + 4 -^ 2 X 1 Xrj 

+16 -^1X1 Xv + 16 -^0X1 Xv ^^61^ — 8/3 (3 x»7 + 2 Xe + X4) -^4X1 

—16/9 (3 Xr; + 2 Xe + X4 + 12 Xi) + 5 X 1 + 64/3 (xe — X 4 ) -^ 7 X 1 -R^ei 

+32/3 (xe + 2 X 4 ) + 6+1 ~ 64/3 (xe + 2 X 4 ) -^ 4 X 1 

-32/9 (Xe + 2X4) LlXi '^ 16)7 ^r]61 + 32/9 (Xe + 2X4) LlXi ^lir) ^r]61 + 16/3 ( 2 X 6 + X 4 ) Llxi) 
+-S'^(Xi) (I 6/9 LgXi Xe-^ i4j? + 32/9 LgXi X 4 -^ler? — 64/3 LgXi-Ri46^ — 32/3 LgXi-^ 145 ,, 
+16/9 1 / 7 X 1 Xe -^i 4 » 7 ^ + 64/9 1 / 7 X 1 X4 Xe Ruri -^ler? + 64/9 + 7 X 1 X 4 R-ieri^ ~ 32/3 + 7 X 1 Xe R-uri 
—64/3 + 7 X 1 X4 Ri6r] + 16 + 7 X 1 + 40/3 -^ 5 X 1 Rueri ~ 3/9 + 5 X 1 Xe Ru-q ~ 16/9 + 5 X 1 X4 Ri^q 
+16/3 + 5 X 1 -Ri 4 e »7 ~ 3 L^Xi Rimq ~ 3 -^ 0 X 1 Rimq ~ 16/3 (xe + 2 X 4 ) -^eXi -Ri 4 e»? 

+16/9 (xe + 2 X 4 ) -^eXi Xe Riai^ + 32/9 (xe + 2 X 4 ) -^eXi X4 Rwq^ 

—32/3 (xe + 2 X4) -^eXi Riisq + 3 (xe + 2 X4) -^4X1 Riieq ~ 16/9 (xe + 2 X4) +4X1 Rwq^ 

—3/9 (xe + 2 X4) +4X1 Riiq^ + 16/3 (xe + 2 X4) +4X1 Riieq) 

+B^ixu Xq) ( - 64/9 Llxi xl R{,q i?^ei + 64/9 L^xi xl Rl,q R^q,, + 64/3 L^xl Rqei^ 

+32/9 Xe Ri4q Rq&i — 32/9 + 5 X 1 X4 Ri&q Rq6i ~ 32/3 -Z 5 X 1 -R^ei 

-64/9 (xe - X 4 ) +/Xi Xe Ruq Rqdi “ 123/9 (xe - X 4 ) +/ Xi X4 ^ler, Rqsi 
+64/3 (xe — X 4 ) + 7 X 1 Rq6i — 64/9 (xe + 2 X 4 ) -^eXi Xe Ruq Rqei 
+64/9 (xe + 2 x 4 ) ^ 6 X 1 X 4 Rl6q Rq6l + 64/3 (xe + 2 X 4 ) L^xl Rq&l 
—32/3 (xe + 2 X 4 ) -^ 4 X 1 Rq&i^ ~ 32/9 (xe + 2 X 4 ) -^ 4 X 1 Ri&q Rqei 
+32/9 (xe + 2 X 4 ) Llxi Rl^q ^^ei) + iXq) (64/9 (xe - X 4 )' ^/Xi Rq,i 
—16/3 (xe + 2 X 4 ) -^ 4 X 1 Xq Rq6i^ + 32/9 (xe + 2 X 4 ) (2 Xe + X 4 ) -^eXi Rqei^ 

— 16/9 (2 Xe + X 4 ) Llxi Xq Rqei^ + 32/9 (2 Xe + xl) Llxi Rq^i^) 

+C''^(Xi) ( — 32/3 LgX? RiiGrj + 16/3 + 5 X 1 -Ri 4 er 7 — 32/3 (xe + 2 X 4 ) -^eXi Riisq 

+16/3 (xe + 2 X4) +4X1 Ri46q) + ^^3(Xi) (3 +3X1 -Ri4e»? ~ 12 -^2X1 — 24 +4X1 + 3 -ZqXi Ri46q) 

+AUxii) ( - 24L5xi - 48LSxi) + A';(Xie) ( - 12+Sxi - 24LSxi) 

+4l2';(X4) ( - 36 L^xi) + A^3(X4e) ( - 48 L^^Xi) + A^Xq) ( - 16 Llxi Rq,i" - 12 

— 16 +0X1 Rq%l ) + -S^(xi) (8 -Z3X1 -Rl 4 e»? + 3 +0X1 -Rl4er7) ( 41 ) 
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= +A{xi) (xi) (xi + 1/9 Xi RIag^^) 

+A{xi)A^{Xu) (4/45 

Rl46r] + 8/9 xi -^146?? + 4/45 (X4-Ilxi) 

+A{xi)A^{Xi,) (2/45 Rl46ri + 4/9xi -Ri46^ + 2/45 (X 6 -Ilxi) -Ri4^^ 

+Z(xi) A^(x46) ( - 2/27x1 Rte,' + 4/27 Xi Rt,, R{,^ - 2/27 Xi R{,,^) 

+^(xi) A^ {Xri) ( - 2/9 Xl Rl46ri R/,61^) + "4(Xl) -^^(Xl) (l/9Xl RuGriRuGr, 

+Xi + 2/9 Xl Ri46ri^) + "4(xi) B^ixi, Xv) ( ~ 8/9 Xi Rugv -^/ei^) 

+Z(xi) C'"(xi) (2/9 xi R{ 4 ,, Rl 4 ,,) + ^(X14) ^^(Xi) (4/45 + 8/9 Xi R\,,, 

+4/45 (x 4 — 11 Xl) Ri&ri) + ^(+ 14 ) A^{X 14 ) ( — 4/27 -R^g^ + 10/27 xi -Rig,, 

+1/27 (2 X, - 2 X4 - 47Xl) + 2/27 (4 x, + 4 X4 + 9 Xi) R/,ei" + 2/27 (4 x, + Xi) R/^ei 
—1/27 (4 X4 + 13 Xl) -Ri 46??) + ^(Xi4) ^^(Xie) ( — Xi) + ^(Xw) ^'^(x»?) ( ~ 2/45 {xrj — X 4 ) 
—8/45 (xv + 9 Xl) -R^ei^ ~ 4/45 (2 x?? — X4 + 9 xi) -R^gi) + ^(Xi4) -R^(xi) (8/9 Xi -Rue?? 

“4/9 (x4 + 2 Xl) Xl -Rie??) + ^(Xu) -R^(Xi? Xv) ( “ 4/9 (x4 + 2 Xi) Xi -R^ei) 

+Z(xi6)7l'"(xi) ( 2/45 -^146r? + 4/9xi -^146?? + 2/45 (xe-llXi) -Ruj?) 

+^(Xi6) ^'^(Xi4) ( — Xl) + ^(Xie) 4l^(xi6) ( “ 2/27 -Ri4g^ + 5/27 xi -Ru?? 

+2/27 (2 X, - 2 xe - 5 Xi) + 1/27 (4 x, + 4 xe + 9 Xi) ^R^gi' - 2/27 (4 x, + Xi) ^^ei 
-1/54 (4 xe + 13 Xi) ^ue??) + ^(Xie) 7l'^(x??) ( - 4/45 (xr? - Xe) 

—4/45 (x?? + 9xi) -R^gi^ + 4/45 (2xr? — Xe + 9xi) -R^ei) + ^(Xie) -R^(Xi) (4/9 xi -Rue?? 

“2/9 (xe + 2 Xl) Xl -Ru??) + ^(Xie) B^iXi: Xv) (4/9 (xe + 2 Xi) Xi -R^ei) 

+^(X4) B^{xi) ( 1/3 Xl X4 -Rie,,^) + ^(X4) -R^(xi? Xv) (2/3 Xi X4 -Rigr? -R^ei) 

+^(X4) B^[xv) ( 1/3 Xl X4 -R^ei^) + ^(X4e) ^^(xi) ( “ 2/27 xi -Rie??^ + 4/27 xi -Rur? -Rie?? 
“2/27 Xl -Ru,,^) + ^(X4e) ^'^(x??) ( “ 2/3 Xi -R^gi^ “ 4/9 Xi -Rig,, -R^gi + 4/9 Xi -Ru)?-R^ei) 
+Z(x4e)i?'"(Xi) ( 4/27 (xe + X4 + Xl) Xl R\4v RlQr) + 2/27(x6“Xi)Xi -Ri4^ 

+2/27 (x 4 “ Xl) Xl -Rie??^) + ^(X4e) B^{xi,Xv) ( “ 4/27 (xe — X 4 “ 3 xi) Xi -Rur? R/^ei 
—4/27 (2 Xe + X 4 + 3 Xi) Xi -Rie?? -^^ei) + ^(X4e) -R'^(Xr?) ( “ 2/9 (3 Xr? + X 4 ) Xi -R^ei^) 
FA{xv) ^^(Xi) ( “ 2/9 Xl -Ruer?-R^ei^) + ^(x??) ^^(Xu) ( — 2/45 (x?? “ X 4 ) 

“8/45 (x,? + 9 Xl) R/^&i - 4/45 (2 x,? “ X4 + 9 xi) -R^gi) + 7l(x,?) 7l'^(xie) ( “ 4/45 (xr? “ Xe) 
“4/45 (x?? + 9xi) -R^gi^ + 4/45 (2xr? “ Xe + 9xi) -R^ei) + ^(x??) ^'^(Xie) ( “ 2/3 Xi -R^gi^ 
“4/9 Xl -Rie?? -Rr?ei + 4/9 Xi -R 14 ?? -^??ei) + ^(x??) A^{xv) (4/9 Xi -R??ei ) 
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+^{Xv) ~ 2/9 Xl Rl46r] + 2/27 Xl Xe Rli-q^ + 1/27 Xl X4 Rwrj^ 

-4/9 Xl Ri46v -^^61^) + MXn) R^ixi, Xn) ( ~ 8/27 Xi Xe Rurj R^ei + 2/27 Xi X4 Rwrj R^ei 
+8/9 Xl + 4/9 Xl -Ri 46»7 + ^(Xr?) R^iXrj) (l/27 (8 xe + X4) Xi R^ei^) 

+Mxv) C^(xi) ( ~ 4/9 Xl Rl46r,R^61^) + -S(Xl) "4^(xi) (1/9 Xl Rl46r,Rl46r, + +1 
+2/9 Xl -Ri46j?^) + Rixi) ^^(Xl4) (8/9 Xl Rl46ri ~ 4/9 (x4 + 2 Xl) Xl -Rl6r?) 

RR{Xi) ^'^(Xie) (4/9 Xl Ri4%ri ~ 2/9 (x6 + 2 Xi) Xi -^ur?) + R{Xi) ^'^(X4) (l/3 Xi X4 Ri%^^ 
RR{Xi) ^^(X46) (4/27 (X6 + X4 + Xl) Xl -^14?? RlGri + 2/27 (X6 ” Xl) Xl Ruri^ 

+2/27 (x4 — Xl) Xl -^16??^) + R{Xi) ^^(Xr?) ( — 2/9 Xl -^140^ -R^ei^ + 2/27Xi Xe 
+ 1/27X1X4-^40^ — 4/9 Xl-R i40^-R ,,01 ) + -B(xi)-S^(xi) (1/9 Xl-^146?) + 4/9 Xl-Rl46r?-^Mer?) 
+Bixi) B^iXuXr,) { - ^/9xiRl4,r,R^,,i") +nxi)c^ixi) (2/9x?i?W) 

+RiXi: Xv) ^^(Xi4) ( — 4/9 (x4 + 2 Xl) Xl R^ei) + RiXi: Xv) ^'^(Xie) (4/9 (xe 
+2 Xl) Xl R^ 6 i) + R{Xi, Xv) ^'^(X4) (2/3 Xl X4 -R10,, Rvei) 

+B{Xi, Xv) ^^(X46) ( - 4/27 (xe - X4 - 3 xi) Xi ^i4r, ^/ei “ 4/27 (2 xe + X4 

+3 Xl) Xl Rwv Rvgi) + B{xi: Xv) ^iXv) ( ~ 8/27 Xi Xe Ruv RvQi + 2/27 Xi X4 -Rie*? -^^ei 

+8/9 xl - 4/9 xl Rl 4 ,v Rv^i) + ^(Xi, X.) ^"^(Xi) ( - 4/9 xl i??4e. ^^ei^) 

+B{Xv) ^^(Xi) (1/3 Xl X4 -R^ei^) + B{xv) ^'^(Xie) ( “ 2/9 (3 x»? + X4) Xi -R^ei^) 

+B{xri) (xri) (1/27 (8 Xe + X4) Xi -R^ei^) + ^'(xi) ^'^(Xi) (2/9 Xi -Ri46»? Rmgv) 

+C{xi) A^(xv) ( - 4/9 xl Rl46v Rv^i) + ^(Xi) B^ {xi) (2/9 xl Rl 4 ,v) 

+^^(xi) ( ~ 1/6X1 Rlfiv + 1/12 Xl -Rl4r?) + "4^(xi)^ (1/2 Xl + 1/18 Xl Rl4fiv^ 

+^^(Xi) A'^{X 14 ) (4 /45 -Ri 40^ + 8/9 Xl Ri4Qv + 4/45 (x4 — H Xi) -^ler?) 

+^^(Xi) 2l'^(xi6) (2 /45 -Ri 40^ + 4/9 Xl Ri 4 Qv + 2/45 (xe — H Xi) -Ri4r?) 

+7l'"(Xi) 71 '"(x 46) ( - 2/27x1 Rl,^^ + 4/27xi R{4v ^le, “ 2/27 xi Rl4v) 

+A^{Xi) A^ixv) ( - 2/9 Xl R‘i4,v Rv^i) + 4l^(Xi) B^ {xi) (l/9 Xi Rl4,v ^ue. + xl 
+2/9 Xl -Ri4ej?^) + ^'^(Xi) -B^(xi, Xv) ( “ 4/9 Xi R^iiv Rv&i') 

+^'^(Xi) C'^(Xi) (2/9 Xl Ruqv Ri 4 &ri) + A^{X 14 ) YBiB ^ ~ -^i46»? + 1/36 (3x77 + 18 Xe 

+42 X4 — Xl) Xl + 1/18 (6 Xr7 + 3 X4 + Xi) Xi -R^ei + 1/18 (6 Xr? + 3 X4 + Xi) Xi -R^ei^ 

+1/18 (3 X4 + 7 Xl) Xl Ri 6 v ~ 1/36 (3 X4 + 7 xi) Xi -^ue??) + A^(xu)"^ ( — 2/27-R140J, 

+5/27 Xl RiQr) + 1/54 (2 Xr7 — 2 X4 — 47 Xi) + 1/27 (4 x?? + 4 X4 + 9 Xi) -R^ei^ 
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+1/27 (4 x-n + Xi) -R^ei ~ 1/^4 (4 X4 + 13 Xi) -^ue??) + ^^(Xu) ^'^(xie) ( — Xi) 

+A^{Xii) A^(Xv) ( - 2/45 (x^ - + 4 ) - 8/45 (x^ + 9 Xi) “ 4/45 (2 Xr, - X4 + 9 Xi) ^/ei) 
+^^(Xi4) -B'^(xi) (8/9 Xi -Ruer? ~ 4/9 (x4 + 2xi) Xi -Rie??) 

+^^(Xi 4 ) -R^(xi) Xv) ( “ 4/9 (x 4 + 2 xi) Xi -R^ei) + ^'^(Xie) ^ 5 /+ ( ~ -Ri 46 ?? 

—1/18 (6 Xr; + 3 xe + Xi) Xi -R^ei + 1/36 (6 Xr? + 3 xe + Xi) Xi -R^ei^ + 1/72 (12 Xr] + 15 Xe 
+36 X4 — Xi) Xi + 1/36 (3 Xe + 7 Xi) Xi -^ 14 ^ — 1/72 (3 xe + 7 Xi) Xi Ruerj) 

+^'^(Xi 6 )^ ( ~ 1/27 -Ri46^ + 5/54 xi -R^4^ + 1/27 {2xri — ‘^X& ~ 5xi) + 1/54 (4x,7 + 4x6 
+9 xi) -R^ei^ “ 1/27 (4 Xv + Xi) -R^ei ~ 1/108 (4 xe + 13 Xi) -^ 140 ^) 

+^'^(Xi6) A^{xr}) ( — 4/45 {xri — Xe) — 4/45 {xv + 9xi) -R^ei^ + 4/45 (2xr) — Xe + 9xi) -R^ei) 
+^'^(Xi6) -R'^(xi) ( 4/9 Xi -Ri 46 ^ — 2/9 (xe + 2 Xi) Xi + ^^(xie) B^ixi, Xv) (4/9 (xe 
+2 Xi) Xi -R^ei) + ^'^(X4) -R'^(xi) ( 1/3 Xi X4 -Rie^^) 

+^'^(X4) B^{xi, Xv) ( 2/3 Xi X4 -R^g^ -R^ei) + ^^(X4) B^{xv) (l/3 Xi X4 -R^ei^) 

+71^(X46) A^ixv) ( - 2/3 Xi - 4/9 Xi T^^ei + 4/9 Xi Rl^v Rv^i) 
+7l^(X46)i?''(Xi)(4/27(x6 + X4 + Xi)Xi RiAv Ri6v + 2/27 (xe - Xi) Xi Ri4v 
A^ jTl (x4 — Xi) Xi -Rie^^) + ^^(X46) -R^(xi, X??) ( ~ 4/27 (xe — X4 — 3 xi) Xi -Ri 4 ^ -R^ei 
—4/27 (2 Xe + X4 + 3 Xi) Xi -Rig^ -R^ei) + ^^(X46) -R^(x»?) ( ~ 2/9 (3 X)? + X 4 ) Xi -R^ei^) 
+^^(Xr?)^ ( 2 / 9 x 1 -R^ei^) + A^{xv) B^{xi) ( — 2/9 xi RimvRv^^ 4” ^/^^^ +1 X6-Ri4^^ 

+1/27 Xi X4 -Rig^^ — 4/9 Xi -Ri4e^ -R^ei^) + ^iXv) B^{xi, Xv) ( ~ 3/27 xi Xe -Ri4^ -R^ei 
+2/27 Xi X4 -Rig^ -R^gi + 8/9 Xi -R^ei^) + A^ {xv) B^{Xv) (l/27 (8 Xe + X 4 ) Xi -R^ei^) 

RA^{Xv) R'^(Xi) ( ~ 4/9 Xi -Ri4g^ -R^gi ) + -B^(xi) ( 1 /I 8 Xi -Ri4e^ + 2/9 Xi -Ri4g^ -R;^ 4 gj^^ 
+B^iXi) B^ixi, Xv) ( - 4/9 xl RIa^v ^+ 1 ') + B^ {xi) ixi) (2/9 xl R{a,v) 

+^^3(Xi) ( 3/4 Xl - 1/2 Xl -Rie^ - l/4xi -R44^) + ^^3(xi4) (l/6xi - 1/3 Xi -R^ei 

—1/3 Xl -R^ei^ ~ 1/3 Xl -Ri6^ + 1 / 6 x 1 -Ri 4 g^) + ^^3 (xi6) ^ 5 ^ (l /12 Xi + 1/3 Xi -R^ei 
—1/6 Xl -R^ei^ ~ 1/6 Xl -Ri4^ + 1/12 Xi -Ri 4 e^) + 77^(1, Xi, Xi) Xi) Xi) (l/3Xi + 2/9 Xi -Ri4g^^) 
+77'^(1, Xl, Xl, Xr,, Xl) ( - 8/9 Xl Ri46v Rv61^) + 77'^(1, Xi, Xu, Xu, Xi) ( - 11/18 Xi 7?u6r, 
-5/6 Xl 72^46^ - 1/9 (4x4-7 xi) Xi 7?i6^) + 77^(1, Xi, Xie, Xie, Xi) ( - 11/36 Xi 7?i46^ 

-5/12 Xl Ri46v - 1/13 (4 Xe - 7Xi) Xi 7?i4^) + 77^(1, Xi, Xr,, Xr,, Xi) (3/9 Xi Rv6i^) 
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+^'^(1, Xi4, Xi4, Xr?, Xi) (l/27 (Xr, “ X 4 ) {Xv “ X4 “ 6 Xi) + 4/27 (xr, + 2 Xi) iXv - 4 Xi) 

+4/27 (xj - X4 X»? - 4 Xi Xr, + Xi X4 - 6 x?) 7?^6i) + Xie, Xie, Xr,, Xi) (2/27 (x,, 

-Xe) iXv - X6 - 6 xi) + 2/27 (x^ + 2 xi) (x»? - 4 Xi) “ 4/27 (xj - Xe Xr? - 4 Xi Xv 
+Xi Xe ~ 6 Xi) R^ei) + 77^(1, X4, Xi4, Xi4, Xi) (3/4 Xi X4) + 77^(1, X46; Xi4, Xie, Xi) (l/2 (xe 
+X4) Xi) + 77^(1, Xvi Xi4, Xi4, Xi) (1/12 Xi Xv + 7/3 Xi Xv T^^ei + 7/3 Xi Xv 77^61^) 

+77'^(7, Xv^ X 16 , X 16 , Xi) ( 7/6 Xi Xv - 7/3 xi Xv 77^61 + 7/6 Xi Xv Rv^i) 

+77'^(2, xi, Xi) Xi) Xi) ( 4/9 Xi 77i4g^ 7?44g^^ + 77^(2, Xi, Xi) Xr?, Xi) ( ~ 8/9 Xi 77i4g^ 77^61^) 
+77'^(2, Xi, Xi4, Xi4, Xi) ( - 5/6 Xi 7744g^) + 77^(2, xi, Xie, Xie, Xi) ( - 5/12 Xi 7744g^) 

+77'^(5, xi, Xi) Xii Xi) ( 2/9 Xi 7744g^^) + (1, Xi, Xi4, Xi4, Xi) (4/9 Xi 7744g^ + 4/3 Xi 7744g^ 

+4/9 (x4 — 4 xi) Xi 77ig^) + 77 /"(1, xi, Xie, Xie, Xi) (2/9 Xi 77i4g^ + 2/3 Xi 7744g^ 

+2/9 (xe - 4 Xi) Xi 77i4^) + (1, Xi4, Xi4, Xr,, Xi) (4/9 {xv “ X 4 ) Xi + 76/9 {xv 

+2 Xi) Xi 77^ei^ + 8/9 (2 x?? — X4 + 2 xi) Xi 77^6i) + Ri (7, Xie, Xie, Xvi Xi) (8/9 {xv ~ Xe) Xi 

+8/9 ixv + 2 Xi) Xi 77/gi^ — 8/9 (2 x,? — Xe + 2 Xi) Xi 77^ei) 

+Hi (2, Xi, Xi4, Xi4, Xi) ( 4/3 Xi 7744g^) + H\ (2, Xi, Xie, Xie, Xi) (2/3 Xi 7744g^) 

+77^(1, Xi, Xi4, Xi4, Xi) (xi 77ig^ — 1/2 Xi 7744g^) + 77^(1, xi, Xie, Xie, Xi) (l/2 Xi 77i4^ 
-l/4x?77^4g^) + 772^1 (l,X4,Xi4,Xi4,Xi) (9/4 x?) + 772^1 (1, Xie, Xu, Xie, Xi) (3x?) 

+77^(1, Xvi Xi4, Xi4, Xi) ( 1/4 Xi + Xi 77^61 + +1 77^61^) + 77]/(1, Xvi Xie, Xie, Xi) (l/2 Xi 
~Xi 77^ei + 7/2 Xi 77^ei^) + 77]/(2, Xi, Xu, Xu, Xi) ( “ 7/2 x/ 7?44g^^ 

+ 77 ]/(2, Xi, Xie, Xie, Xi) ( - 7/4 Xi 77/4g^) + 772^(1, Xi, Xu, Xu, Xi) ( - Xi 77/g^ 

+7/2 Xi 7744g^) + 77]/(l, xi, Xie, Xie, Xi) ( ~ 7/2 xi 7?44^ + l/4xi 7744g^) 

+772^(7,X4, Xu, Xu, Xi) (-9/4xi) + 772^(7, Xie, Xu,Xie, Xi) (-3Xi) 

+772^(7, Xvi Xu, Xu, Xi) ( - 7/4 xi - Xi 77^ei “ Xi 77^ei^) 

+77]/(l, Xvi X 16 , Xie, Xi) ( ~ 7/2 xi + Xi 77^61 ~ 7/2 xi 77^61^) 

+77]/(2, xi, Xu, Xu, Xi) ( 7/2 Xi 77/4g^) + 77]/(2, xi, Xie, Xie, Xi) (7/4xi 77/4g^) (42) 


B Expressions for the decay constant 

= +A^{xia) ( 7 ) + A'"(xie) (7/2 


(43) 




+A^(Xi) ( - 4/3 LI XI FU, + 4 1, + 4 Xi R^,,, - 70 Xi - 4 L/ Xi 
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+4 1, + 4 LI XI Ruev) + (4 LI - 4 (xe + 2 X4) t, - 10 (x4 + Xi) ^3 

—4 (x4 + Xi) -^ 0 ) + ^^(Xie) ( 21/5 Xi “ 2 (xe + 2 X 4 ) L^ — 5 (xe + Xi) Ll — 2 (xe + Xi) -^ 0 ) 
+A'"(X4) (12 LIxa-Q L; X4 - 24 LI X 4 ) + A'"(X46) (8 (xe + X 4 ) L^ - 4 (xe + X 4 ) L; 

-16 (X6 + X4) L[) + A^ix,) (8/3 LI xi x, R^" - 2 x, - 8 x. 


-8 LI X, RW + 4/3 (2 X6 + X4) LA + B^Xi) ( - 4/3 L; xi R^,, + 4 xi R^e, 


+4 Lq Xi -Ri 46??) + -S^(Xi4) (8 (x4 + Xi) (Xe + 2 X4) L^ — 4 (x4 + Xi) (Xe + 2 X4) 

+4 (X4 + Xi)^ LI — 2 (x 4 + Xi)^ Z's) + R^iXw) (4 (X6 + Xi) (X6 + 2 X4) /'e 

“2 (xe + Xi) (Xe + 2 X4) + 2 (xe + Xi)^ L^ — (xq + Xi)^ Z's) + ^^3(Xi) ( ~ 41/3 -^146*? 


+6 LI+ 12 LI-A LI + A^^Xu) (12 + 24 LA + A^.^Xw) (6 L; + 12 L 


tV' 


iV" 


+^23(X4) (181/2] + ^23(X46) (24//^ j + A 23VA»? 
’ 23 IX 1 ) (“41/3 -Ri46^ — ALq -Ri4g^ 


8L5/R^642 + 6L^2 + 8LSi?^6i^ 


+5^ 


(44) 


^42 ^ (l/18/Rf46. - 1/9 (X4 + 2xi) R^r,) 

+A{xi)B^{Xie) (l/36 Rl46ri 1/18 (x 6 + 2 xi) -^14r?) 

+Z(xi4) /4 '"(xi4) ( 5/54 - 5/27/R^gi - 3/27R^,^ - 3/27R^^ + 5/54/?^,) 

+A(xi 6 ) /4'"(xi 6) (5/108 + 5/27i?^64 - 5/54 _ 5/54^.^^ + 5/108 72^,) 

+A{xA B^Ixm) ( — 1/36 (xr? — X4) — 1/9 [xt] + X4 + Xi) R/^ei ~ 1/9 iXv ~ Xi) -R^ei^) 
+A{xA) ixio) ( ~ 1/18 (Xr? — Xe) + 1/9 (Xr? + Xe + Xi) -R^ei “ 1/18 {Xv ~ Xi) -R^ei^) 
+-B(X14) A'^ {xi) (I/I8 -R44g^ — 1/9 (X4 + 2 Xl) Rl&ri) + -S(Xl4) {Xr/) ( ~ 1/36 (Xr? “ X4) 

~l/9 (Xr; + X4 + Xl) -R^ei ~ 1/9 (Xr; ~ Xl) -^^61^) + -®(Xie) /4^(xi) (l/36 -Ri4g^ 

~1/18 (xe + 2 Xl) -Ri4r?) + -S(Xie) ^{Xt/) ( — 1/18 (xr? — Xe) + 1/9 iXv + Xe + Xi) R/jei 
-1/18 (x, - Xl) + 4l'"(Xi) (l/8 Xl - 1/12 Xl RAv “ 1/24 Xi RA^) 

+^'^(Xi) R^iXu) (l/18 -Ri4g^ — 1/9 (x4 + 2 Xl) -Rierj) + ^^(Xi) -S^(xie) (l/36 -Ri4g^ 
-1/18 (xe + 2xi) Rur] Rimr] 1/24 (Xrj + 6 Xe + 14 X4 + 5 Xl) 

—1/12 (2 x») + X4 + Xl) -R^ei “ 1/12 (2 X»? + X4 + Xi) -R^ei^ ~ 1/12 (X4 + 3 xi) -Rfg^ 

+1/24 (X4 + 3xi) + ^""(Xu)' (5/108 - 5/54 ii^g^ - 5/54 - 5/54 i?fg^ 

+5/108 i?44g^] + A^(xig) (l/24-Ri4g^ + 1/12 {2xn + Xe + Xi) -R^ei 

—1/24 (2 Xr; + Xe + Xl) R/^ei^ ~ 1/48 (4 + 5 xe + 12 X4 + 5 xi) — 1/24 (xg + 3 xi) Rur) 
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+1/48 (X 6 + 3xi) (5/216 + 5/54 - 5/108 - 5/108 

+5/216 -Ri 46 ^) + {xr)) ( ~ 1/36 (x^ — X 4 ) — 1/9 (x»? + +4 + Xi) -R^ei 

~l/9 (Xr; — Xi) -R^ei^) + ^{Xr}) {xw) ( — 1/18 (xv ~ Xe) + 1/9 (x»? + Xe + Xi) -R^ei 
-1/18 (x»? - Xi) -R^ei^) + ^^ 3 (Xi) Yq^ ( “ + l/4-Ri6r) + 1 / 8 -^Ur?) 

+^^(Xl4) ~ 1/12 + 1/6 RrjQi + 1/6 R^qi + 1/6 -Rigj) — 1/12 -Ruer?) 

+^^(Xi6) Yg//2 ( ~ 1/24 — 1/6 R^qi + 1/12 R^qi + 1/12 -Ri 4^ — 1/24 -Ri46^) 

+^'^(1, Xi, Xi4, Xi4, Xi) (1/12 Bi4Qn 1/6x1 B-lQri + 1/12x1 -^146r?) 

+i7'^(l, Xi, X 16 , X 16 , Xi) (l/24-R44g^ — 1/12 xi -Ri4^ + 1 / 24 x 1 -^ 446 , 7 ) 

+H^{1, X4, Xi4, Xi4, Xi) ( — 3/8 X 4 ) + -^'^(1, X46, Xi4, X 16 , Xi) ( — 1/4 (xe + X 4 )) 
+^'^(l,X»?,Xi4,Xi4,Xi) ( - 1/24x7? - l/6x,?i?^6i “ l/6Xr?^^6i^) 

+^'^(1,X7?,X16,X16,Xi) ( - 1/12X7? + l/6Xr?^^61 - 1/12X7? ^^61^) 

+B^(2, Xi, Xu, Xu, Xi) (5/54 Xi -Rugt?) + ^^(2, Xi, X 16 , X 16 , Xi) (5/108 Xi i?i 467 ?) 

+BI' (2, xi, Xi4, Xi4, Xi) ( - 1/108 xi ^u 67?) + ^1"(2, Xi, X16, X16, Xi) ( - 1/216 Xi ^mgt?) 

+^ 1 ^^(3, Xi4, Xi, Xi4, Xi) ( - 1/54 xi ^ 1467 ?) + ^ 1 "(3, X 16 , Xi, X 16 , Xi) ( - 1/108 Xi ^U 67 ?) 
+-^^(17 XI 7 Xl4, Xl4, Xl) (1/2 -Ri6,7 ~ 1/4-^1467?) + -^^( 17 X 17 X 167 X 167 X 1 ) (1/4-Ri4,7 — 1/8 

+-^2^(17 X4, Xl4, Xl4, Xl) (9/8) + Rr2^(l, X46, Xl4, XI6, Xl) (3/2) 

+-^^(17 Xr/, Xu, Xl4, Xl) ( 1/8 + 1/2 i?^g 4 + 1/2 i?^g 4 ^) 

+-^^(17 X7?7 XI67 XI67 Xl) (1/4 — l/ 2 i ?^64 + l/ 4 i?^ 64 ^) 

+77^(2, X 17 Xi 4 , Xi4, Xl) ( — 1 / 4 -^ 146 , 7 ) 4” -^^(2, X 17 X 16 , X 16 , Xl) ( — 1 / 8 -^ 446 , 7 ) 

-hff '^(l, Xl, X 17 X 17 Xl) ( 1/6 Xl + 1/9 Xl -Ri467?^) + ^(17 X 17 X 17 X7?7 Xl) ( — 4/9 Xl -R 1467 ? -R^6i^) 

+-^'^^(17 X 17 X14, Xi4, Xl) ( - 11/36 Xl -Ri 467 ? - 5/12 Xl -RUti “ 1/18 (4 X 4 - 7 Xi) Xi -Rig,,) 
+-R'^(l7Xi7Xi67Xi67Xi) ( - 11/72 Xi-Ri46,7 - 5 /24 X? -RUti “ 1/36 (4x6 - 7Xi) Xi -R^,,) 

+-R'^(17 X17 Xt?, X7?7 Xl) (4/9 Xl -R^e/) + 77'^(1, Xi4, Xu, Xr?, Xi) (l/54 (x,? - X4) (Xr? - X4 
-6 Xl) + 2/27 (x,? + 2 Xl) (Xr? - 4 Xi) -R/gi^ + 2/27 (xj - X4 Xr? - 4 Xi Xt? + Xi X4 - 6 Xi) -R^ei) 
+-R'^(1, X16, X16 , Xt?, Xl) ( 1/27 (xt? - X6) (Xr? - X6 - 6 Xl) + 1/27 (x,? + 2 Xl) (Xr? - 4 Xl) 

B^6i^ - 2/27 (xj - X6 Xt? - 4 Xl Xt? + Xl X6 - 6 Xi) -R^ei) + 77'^(1, X4, Xm, Xi4, Xi) (3/8 Xi X4) 
+-R''^(1, X46, Xl4, XI6, Xl) (1/4 (X6 + X4) Xl) + -R'^(1, Xr?, Xl4, Xl4, Xl) (l/24 Xl Xr? 
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+1/6 Xi XrjR'r^ei + 1/6X1 XnKiiii) + ^'^(1, Xr?, X16, X16, Xi) (1/12 Xix»? - 1/6 Xi X^KiGI 
+1/12 Xi Xf) + R ^(2) Xi) Xi) Xi) Xi) (2/9 Xi -Ri46»? -^14677) 

+i7 ^{2, xi, Xi) X»?7 Xi) ( “ 4/9 Xi -^1467? '^(2, Xi, Xi4, Xi4, Xi) ( ~ 5/12 xl RiAGri) 

+H ^(2, Xl, X 16 , X 16 , Xl) ( ~ 5/24Xl -R 14677 ) + H ^^(5, Xi, Xi, Xi, Xi) (l/9 Xi -^ 1467 ?^) 

+i7i^(l, Xl, Xi4, Xi4, Xl) ( 2/9 Xl -^1467? + 2/3 Xl -^ 1467 ? + 2/9 (x4 — 4xi) Xi -^igt?) 

Xl, X16, X16, Xl) (1/9 Xl -^1467? + 1/3 Xl -^14677 + 1/9 (X6 — 4xi) Xl -^1477) 
xi4, Xi4, X77, Xl) (2/9 (xt? — X4) Xl + 8/9 (X77 + 2 Xl) Xi R^ei^ + 4/9 (2 — X4 

+2 Xl) Xl R^ 6 i) + -^1^(17 X16, X16, X77, Xl) (4/9 {xri — Xe) Xl + 4/9 (X77 + 2 Xi) Xi -R/ei^ 

—4/9 (2 X77 — Xe + 2 Xi) Xi R^ei) + -f^i'^(2, Xi, Xi4, Xi4, Xi) (2/3 Xi -^ 14577 ) 

+H^{2, Xl, XI6, XI6, Xl) (1/3 Xl -R14677) + -^2^(1; Xl, Xl4, Xl4, Xl) (1/2 Xl Rwrj “ l/4Xl -R14677) 
+^2^(17 Xi,Xi6,Xi6,Xi) (i/ 4 x?^M 77 “ 1/8 x?^u 677) +^^2^(17X47X147X147X1) (9/8 X?) 

+-^21"(I7 X467 Xi 47 X167 Xl) (3/2 Xl) + (1, Xt ), Xi 47 Xi 47 Xl) (1/8 Xl + 1/2 Xl R^ei 

+1/2 Xl -R^ei^) + -^2r(l7 X777 X167 X167 Xl) (1/4 Xl ~ 1/2 Xl R ^ 6 i + 1/4 Xi R ^ gi ^) 
+^2^(27X17X147X147X1) ( - 1/4 X? ^14677) +^21/(27X17X167X167X1) ( - 1/8 x?^m 677) 

+-^2^(17X17X147X147X1) ( “ 1 / 2 X 1-^1677 + l/ 4 xi -^14677) 

+-^/7(l7Xl7Xl67Xl67Xl) ( - 1/4X1-^1477 + 1/8X1-RU677) +-f^27(l7X4,Xl4,Xl4,Xl) ( “ 9/8 Xl) 
+-f^/7(l7 X467 Xl47Xl67Xl) (-3/2Xl) +i//7(l,X77 7 Xl47Xl47Xl) (-1/8X1 - 1/2 Xi-R /61 
—1/2 Xl R^6 i^) + -^2^(17 Xrj, X167 X167 Xl) ( — l/4xi + 1/2 Xl R^ei ~ 1/4 Xi R^ei^) 
+H2^{‘^,Xi,Xu,Xii,Xi) (1/4X1-R14677) +-^27(27X17X167X167X1) (1/8X1-R14677) (45) 
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